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The quantum theory of antiferromagnetism in metals is necessary for our un- 
derstanding of numerous intermetallic compounds of widespread interest. In 
these systems, a quantum critical point emerges as external parameters (such 
as chemical doping) are varied. Because of the strong coupling nature of this 
critical point, and the "sign problem" plaguing numerical quantum Monte 
Carlo (QMC) methods, its theoretical understanding is still incomplete. Here, 
we show that the universal low-energy theory for the onset of antiferromag- 
netism in a metal can be realized in lattice models, which are free from the 
sign problem and hence can be simulated efficiently with QMC. Our simula- 
tions show Fermi surface reconstruction and unconventional spin-singlet su- 
perconductivity across the critical point. 
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The presence of an antiferromagnetic transition in a metal is common to compounds such as 
electron-doped cuprates (i), iron based superconductors (2), and heavy fermion Kondo lattice 
systems (3). Whereas our understanding of quantum antiferromagnetism in insulators has seen 
remarkable advances (4), analogous problems in metals are far more complicated because of 
the subtle interplay between the low energy fermionic quasiparticles on the Fermi surface, and 
the quantum fluctuations of the antiferromagnetic order parameter. In addition, the presence 
of the Fermi surface has hampered large scale numerical studies, because QMC algorithms 
are afflicted by the well-known fermion sign problem. Such algorithms express the partition 
function as a sum over Feynman histories, and the sign problem arises when the weights as- 
signed to the trajectories are not all positive because of quantum interference effects. A general 
solution to the fermion sign problem has been proved to be in the computational complexity 
class of nondeterministic polynomial (NP) hard (5), and so there has been little hope that the 
antiferromagnetic quantum critical point could be elucidated by computational studies. 

Application of the methods of quantum field theory and the renormalization group to the 
onset of antiferromagnetism in a metal (<5), has identified (7, 8) a universal quantum field the- 
ory which captures all the singular low energy quantum fluctuations that control the quantum 
critical point and deviations from the Fermi liquid physics of traditional metals. In two spa- 
tial dimensions, the field theory is expressed in terms of fermionic excitations in the vicinity 
of a finite number of 'hot spots' on the Fermi surface, and is thus independent of the details 
of the fermionic band-structure, except for the number of hot-spots and Fermi-velocities at the 
hot-spots (9). Recent work (10, 11) has shown that the renormalization group and Feynman 
graph expansions of the field theory flow to strong coupling, making further analytical progress 
difficult. 

Here, we show that the universal quantum field theory can be realized in lattice models 
which are free of the sign problem, and so is amenable to large scale QMC studies. Our claim 
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does not contradict the no-go theorem of Ref. (5), because we do not provide a general recipe 
for eliminating the sign problem. However, we will eliminate it for the specific case of the onset 
of antiferromagnetic order in a two-dimensional metal, provided the perturbative arguments on 
the importance of the hot spots to the quantum field theory (7,8, 10, 11) apply. Our modified 
lattice model has at least two bands. Therefore, in cases in which there is only a single active 
band at the transition, such as in the electron-doped cuprates, our method requires modifying 
Fermi surface far away from the hot spots; we show that this can be done while preserving the 
universal low-energy properties of the antiferromagnetic critical point. On the other hand, our 
method applies to multi-band situations (such as in the iron-based superconductors) without 
changes to their Fermi surface configuration. Being a low-energy effective theory, the method 
will not apply where the proximity of a Mott insulator is important, as is likely the case in the 
hole-doped cuprates (12-16). 

To illustrate our method, we now consider the onset of antiferromagnetic order in a simple 
one-band model on the square lattice, as is appropriate for the electron-doped cuprates. The 
electrons, Ck (the spin index is left implicit), with dispersion e^, have a single "large" Fermi 
surface (Fig UK)- The antiferromagnetic order parameter is y? q ; we will assume the important 
fluctuations of <^ q are restricted to small values of |q|, much smaller than the size of the Brillouin 
zone. The antiferromagnetic ordering wavevector is K = (7r, h), and y? q represents the electron 
spin density at the wavevector K + q; we will also refer to the antiferromagnetic order as spin 
density wave (SDW) order. We can now write the electron part of the Hamiltonian as 



where A is the 'Yukawa' coupling between the electrons and the SDW order, and s are the Pauli 
matrices. The Yukawa term is the simplest coupling consistent with translational symmetry and 
spin-rotation invariance, and can be derived e.g. by decoupling of the repulsive interaction in a 




(1) 



k 



k,q 



3 



Hubbard model by an auxiliary field which maps to (p in the long- wavelength limit (17). The 
hot spots are at k for which e k = £k+K = (Fig. [IK); at these points, v? q= o scatters electrons 
between initial and final states which are both on the Fermi surface. To obtain the electron 
Fermi surface in a metal with SDW order, we replace </? q by its expectation value (<p q ) = N5^ 
(where N is the staggered magnetization), and recompute the electron dispersion; this leads to 
the Fermi surface reconstruction shown in Fig. [I}3. 

We now describe our method to deform the model, such that the sign problem is avoided, 
while preserving the structure of the hot spots. Let us separate the hot spots into two groups, 
so that K only connects hot spots from one group to the other. Now deform the one-band 
electronic dispersion to a two-band model with an additional 'orbital' label, so that all the 
hot spots in one group are on the Fermi surfaces of the first band, while the hot spots of the 
other group reside on the Fermi surfaces of the second band (an example of such a dispersion is 
shown in Fig.QJJ, in which the 'horizontal' and 'vertical' Fermi surfaces are part of two separate 
electronic bands). Note that the vicinities of the hot spots in the two-band model are essentially 
identical to those in the one-band model in Fig. [TK, and so the same low energy theory for the 
onset of antiferromagnetism applies to both models. With no further assumptions, the deformed 
model has only positive weights in a suitable quantum Monte Carlo realization. 

We will write down a specific lattice model for which we will establish a sign-free Monte 
Carlo algorithm, and then present numerical results. We begin with the band structure of the 
c k electrons in Fig. Qjj. We write the band with vertical Fermi surfaces in terms of fermions 
ip x with Ck — > 4>x,\n, and the band with horizontal Fermi surfaces in terms of fermions ip y with 
c k — > ipy,k+K' This leads to the ip x :V Fermi surfaces shown in Fig. |2K- Then our model has the 
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action S = Sp + S v = J Q dr(L F + L v ) with 



i,j,a=x,y i 




Here i,j run over the sites of the square lattice, r is the imaginary time and (3 - the inverse 
temperature. The parameter r will be used to tune across the quantum critical point, and u 
is a non-linear self-coupling of (p. The ip x (ip y ) fermion hops along the horizontal (vertical) 
direction with an amplitude tn = —1 (+1), and along the vertical (horizontal) direction with an 
amplitude t± = —0.5 (0.5), respectively; the resulting band structure is shown in Fig.[2K (solid 
lines). The model has C 4 symmetry, and its apparent violation is an artifact of the shifting of 
the ip y fermions by K. We chose the chemical potential [L\ — [ii — —0.5, c = 1, u = 1, and 
A = 1. 

By construction, the modified two-band model has the same hot spot structure as the original 
one-band model. Therefore, we argue that it preserves the universal properties of the antiferro- 
magnetic transition. We prove (9) that the introduction of the second band eliminates the sign 
problem in this model. 

Note that it is possibly to analytically integrate out (p in Eq. (IS 121) . and establish equivalence 
to a large class of Hubbard-like models to which our method applies. However, we choose 
to keep (p as in independent degree of freedom because it keeps the physics transparent and 
streamlines the analysis. 

We have performed determinant Monte Carlo simulations of the action (IS 121) using the algo- 
rithm described in Refs. (18-20), for systems of linear size up to L = 14 and inverse temperature 
(3 = 14, with either periodic or anti-periodic boundary conditions. An imaginary time step of 
At = 0.1 was used in most of the calculations; we checked that the results do not change for 
At = 0.05. Up to 50000 Monte Carlo sweeps were performed for each run, giving a statistical 
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error for most measured quantities of a few percent. 

First, we present results showing the reconstruction of the Fermi surface across the SDW 
transition. Fig. |3] shows the fermion occupation number summed over the two flavors of 
fermions as a function of quasi-momentum. The Fermi surfaces are clearly visible as disconti- 
nuities, r = 0.5 is found to be on the disordered side of the SDW critical point, and the Fermi 
surface closely resembles the one in Fig. |2]A At r = 0, a gap opens at the hot spots, and the 
Fermi surface is reconstructed into electron and hole pockets, as in the SDW ordered state in 
Fig. [2b. Decreasing r further to —0.5 increases the magnitude of the SDW order parameter, and 
causes the hole pockets to disappear and the electron pockets to shrink. 

To examine the magnetic transition, we computed the SDW susceptibility x^ — Si jo dr (tpi (r ) • 
y? (0)). Figure SK shows x<p normalized by L 2 /3 as a function of r. In order to extract informa- 
tion about the zero-temperature limit, we scale (3 with the linear system size; in the appropriate 
units, (3 = L was used. We observe a rapid upturn in Xp near r = 0.25. For r < 0.25, x<pl (L 2 f3) 
for different system sizes and inverse temperatures nearly collapse on top of each other, which 
is the expected behavior on the ordered side of the transition. The results are consistent with 
a second-order transition at r c « 0.25. This is further supported by the Binder cumulant in 
Fig. |4}3, where we observe the expected behavior in both phases, separated by a critical point 
atr c = 0.25 ±0.1. 

The SDW critical modes mediate effective inter-fermion interactions, which can lead to in- 
stabilities of the Fermi surface. As a result, additional competing phases can appear. Near 
the SDW critical point, these instabilities are a result of a subtle competition between the en- 
hancement of the SDW fluctuations, which tends to strengthen the effective interactions, and 
the loss of coherence of the fermionic quasi-particles {10, 11). Superconductivity is a natural 
candidate for the leading potential instability. In order to examine the emergence of a supercon- 
ducting phase near the SDW critical point, we have computed equal-time pairing correlations 
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P±{xi) = (A ± (f i )A ± (0) t ). Here, A ± (£;) = is v ah {^ ixa ^ ixh ± ipiy a ipiyb) (where a,b =t4 are 
spin indices) are superconducting order parameters with either a + or — relative sign between 
the two fermionic flavors (square lattice symmetry A lg and B lg , respectively). 

In order to probe for long-range order, we measured P±(xi) near the maximum range 
Xmax = (L/2,L/2). We plot P±{x rnax ) = ^Y, ex , y =-i P ±(x™<™ + e ^Vx + e yVy)> where ff x = 
(1, 0) and fj y = (0, 1), in Fig. H£. Long-range superconducting order at /3 — > oo would corre- 
spond to superconducting correlations that saturate to a constant upon increasing L and f3. 

The Bi g pairing correlations are found to be significantly enhanced in the vicinity of the 
SDW critical point, r c w 0.25. The A lg correlations are significantly smaller in magnitude and 
negative in sign. This is consistent with the expectation that the effective attraction mediated by 
magnetic fluctuations promotes superconductivity with a sign change between the two orbitals 
(22,23). 

The maximum of the B\ g correlations occurs for r ps 0.5, on the disordered side of the 
magnetic critical point which is located at r c w 0.25 (21). Interestingly, the suppression of the 
superconducting correlations away from the optimal r is very asymmetric: whereas the pairing 
correlations decrease gradually for r > r c , they are suppressed dramatically for r < r c . This 
may be a result of the opening of an SDW gap on portions of the Fermi surface. 

The method described in this Letter opens the way to study various physical aspects of 
spin density wave transitions in metals, in a numerically exact way. The interplay between 
unconventional superconductivity and magnetism and possible non-Fermi liquid behavior in 
the quantum critical regime should now be accessible. Moreover, such simulations will provide 
controlled benchmarks for analytic approximations (7, 8, 10, 11). 

The two-band model presented here is a member of a wider family of strongly correlated 
fermionic models that can be rendered free of the sign problem. It has already been known that 
some models with two flavors of fermions interacting via a four-fermion interaction are sign 
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problem free at generic fermion density (24). Remarkably, these models do not rely on any 
specific characteristic of the electron dispersion; e.g. there is no requirement for particle-hole 
symmetry, or any symmetry that relates the two bands. Extensions of this trick to related models 
of physical interest should be possible. 
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Figure 1 : (A) Fermi surface of the Fermi liquid phase of a single band model on the square 
lattice with unit lattice spacing. The "hot spots" are denoted by the filled circles. (B) The 
reconstructed Fermi surface in the metal with SDW order. The dashed lines show the Fermi 
surface in the metal without SDW order, and its translation by K. Gaps have opened at the 
hot spots, leading to small "pocket" Fermi surfaces. (C) A deformed Fermi surface of the 
metal without SDW order, in which the vicinities of the hot spots are unchanged from (A). The 
horizontal and vertical Fermi surfaces now belong to separate electronic bands. 
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Figure 2: (A) Fermi surfaces (full lines) of Lp for free ip XiV fermions with the parameters listed 
in the text. The dashed lines show the portion of the Fermi surface in Fig. [TJ; which was shifted 
by K to obtain the ip y Fermi surface. Now the hot spots are at the intersections of the Fermi 
surfaces. (B) Mean-field ip XtV Fermi surfaces with SDW order | (</?) | = 0.25. 
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Figure 3: Quantum Monte-Carlo results for the fermion occupation number = (ipl k ip x k + 
VykVVk) /2 as a function of k across the Brillouin zone, for systems with L — 14, /3 = 14, 
and r = —0.5, 0, 0.5. In order to enhance the resolution, results from simulations with either 
periodic or anti-periodic boundary conditions in the x and y directions were combined. Despite 
appearances, full square lattice symmetry is preserved in all our computations for the original 
Ck fermions. 




Figure 4: (A) The SDW susceptibility x<f> normalized by L 2 f3, as a function of r, for systems 
of size L — 8, 10, 12, 14 and (3 = L for each curve. The statistical errors in \<t> are smaller than 

the symbol size. (B) The Binder cumulant for an 0(3) order parameter Cb = 1 — ^~pp> where 

$ = -h J2i 0%i approaching the expected values of 0.4 and in the two phases. (C) Equal-time 
pairing correlations in systems of size L = 10, 12, 14 and (3 = L for each curve, as a function 
of r. Dashed (solid) lines show P + (P_), corresponding to A\ g {B\ g ) superconducting order 
parameters, in which the pairing amplitude in the two fermion flavors is of the same (opposite) 
sign, respectively. r c is the estimated position of the SDW critical point. 
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Supplementary Material 
Quantum Field Theory 

In this section, we briefly review the universal field theory for the SDW transition developed 
in Refs. (7,8,10,1 1). This effective field theory has been analyzed perturbatively in the order 
parameter-fermion coupling A. 

One of the key results of the analysis is that, if we ignore the possibility of a high-energy 
pairing instability, the fermion and order parameter propagators acquire universal singular parts 
which depend only on the structure of the hot spots (e.g. the velocities at the hot spots and the 
angle between them). This justifies the assumptions behind the construction of the lattice model 
presented in the main text: as long as the structure of the hot spots is preserved, we expect the 
universal behavior near the antiferromagnetic critical point to be unchanged. The microscopic 
parameters of the model should only come in through the ultraviolet cutoffs to the critical fluc- 
tuations. These cutoffs can, in principle, be set by matching at high energy scales. Specifically, 
we can match a sign-problem-free lattice model to a Hubbard-like model by equating their hot- 
spot Fermi surfaces and Fermi velocities. Other parameters of the sign-problem-free model can 
be determined by matching its physical observables to those of the sign-problem-present Hub- 
bard model at temperatures high enough to allow accurate computations by other methods for 
the latter model. The sign-problem-free model can then be used to compute observables at low 
temperatures. 

Another observation made in Ref. (10), however, is that ultimately the conventional ways to 
control the perturbative series for the effective field theory, such as an expansion in the inverse 
number of fermion flavors, are uncontrolled for this problem. The fate of the flow to strong 
coupling has to be resolved by numerical simulations. For this we need a lattice regularization 
of the continuum quantum field theory, and the lattice model considered in the main body of the 
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paper provides precisely such a regularization. 

The field theory is formulated in terms of the fermion excitations in the vicinity of the hot 
spots. The antiferromagnetic order parameter <p, with wavevector K, connects fermions at a hot 
spot at wavevector k with fermions at a hot spot with wavevector k + K; both fermions are on 
the Fermi surface if = £k+K = 0, and this defines the allowed values of k. Linearizing the 
fermion dispersion about the hot spots, and expanding the order parameter in spatial gradients, 
in two spatial dimensions we obtain the Lagrangian L = + L v where 

£>i> = ip\ (d T - ivi • V) ipi + i)\ (d T - iv 2 -V)ip 2 + X0- Ol#2 + H.c.) 

= + \mr + \(f + \m (SD 

Here ip a , with a = 1, 2, are two species of low energy fermions in the vicinity of the hot spots 
at k and k + K, and v a are their Fermi velocities. A similar Lagrangian applies to the other hot 
spots. This theory has the same general structure as coupled fermion-boson theory in particle 
physics, such as the Gross-Neveu model (25), with fermions and bosons coupled via trilinear 
"Yukawa" coupling A. The key difference is in the fermion dispersion, which does not have a 
relativistic form. In the relativistic cases, the fermion dispersion has a Dirac form with energy 
~ ±i>|k|, and this vanishes only at isolated points in the Brilluoin zone. The resulting fermion- 
boson theory is well understood (25). In our case, the fermions dispersion ~ v ■ k, and this 
vanishes on a line in the Brillouin zone which is orthogonal to v. This is the central difference 
which makes the quantum field theory in Eq. dSTT ) strongly coupled. 
Let us parametrize the Fermi velocities by 

Vi = (v x , v y ) , v 2 = (-v x , v y ) . (S2) 

Here, for convenience, we have rotated the coordinates by 45° relative to Fig. la in the main 
text. It is useful to introduce the ratio and the modulus 

Vy 

tan0 = — , v = v . (S3) 

v x 
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Here < 20 < ix is the angle between the Fermi surfaces at the hot spot. We now present the 
singular terms in the low energy spectrum at two loop order, as obtained in Refs. (10,17). We 
will restrict the expressions to precisely at the quantum critical point at zero temperature; all 
our conclusions, and similar but lengthier expressions, apply also close to the quantum critical 
point and at low temperatures. For the fermion Green's function we have 



\ . 3usin20. / / (va-p) 2 |v B • p| , 
G a (w, p) = -v a ■ p + 1 sgn(w) I y 7|w| + — - 2 — ) , (S4) 



where to is an imaginary frequency, 1 = 2 and 2 = 1, and 

7 = 7; > (S5) 

where Nh is the number of pairs of hot spots (Nh = 4 for the electron-doped cuprates). In the 
expression (|S4l we have dropped the bare free fermion iui term because it is not as singular as 
the self-energy correction from the fluctuations of the anti ferromagnetic order, and we have not 
explicitly written the real part of the self energy which renormalizes the velocities v x and v y . 
The singular part of the propagator of the boson (p is 

£) _1 (o;,p) = 7|w| + p 2 . (S6) 

In these expressions above the strength of the interactions is controlled by the Yukawa cou- 
pling A, and hence via the value of 7. However a key observation is that dependence on A can 
be scaled away, and the above low energy spectra are actually universal. Indeed, it is easily 
seen from Eqs. (IS4IS6I ) that after rescaling momenta by p — > Ap, the A dependence appears as 
overall prefactors which can be absorbed into a rescaling of the fields. This independence on 
the value of A is a general feature of the low-energy quantum field theory (10), and is crucial 
to its properties. One of its consequences appeared in the leading log estimate of the pair- 
ing instability presented in Ref. (10), which was found to be a logarithm- squared term with a 
A-independent prefactor. 
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Quantum Monte Carlo 



We set up S F + S<p for a Monte Carlo study (17). Discretizing imaginary time, the partition 
function becomes 



dipexp {-S v ) Tr^ 



N 



IP- 

n=l 



O (At 2 ) 



(S7) 



where Ar is an imaginary time step, (3 = NAr, and the operators B n are given by 



B r 



K and V n are matrices given by 



(S8) 



^i,j,a,a' ;s,s' ^s,s'^a,a' ( ^a,ij A*) 
K;i,i;a,a';s, S ' = A 5ij [s • Ifii (nAT)] ss , , (S9) 

where i, j are spatial indices, o\ is a Pauli matrix, a, a' = x,y are flavor indices, and s, s' =|, I 
are spin indices, ip^ is a vector of fermionic operators, 

^ = (V£i jt > • ■ • > ^,1,4.' • • • ' ^,4-' 

^,i,t> • • • ' ^Lv\t' ^J,U> • • • ' ^Lvu) • (S10) 
J\f is the number of lattice sites. Note that Tr^, in Eq. IS7l represents a trace over fermionic 
many-body states in Fock space. This fermionic trace can be carried out, giving 



Tr 



A' 



IP 



_n=l 

At K ArV n „ — ^ At K 



det 



iV 



1 + ; \B r 



n=l 



(Sll) 



where B n = e * e ^ TVn e 2^ TK . For a proof of this formula, see, e.g., Ref. (19). We then 
arrive at the following form of the partition function: 

N 

Z = I dipexp (— S v ) det 1 + Y\ B r , 



n=l 



+ O (Ar 2 ) 



(S12) 



which can be evaluated using Monte Carlo techniques, by sampling over space-time configura- 
tions of (fii (r). 
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Positivity of the action 

Monte Carlo sampling can be done efficiently if the action in Eq. IS12l is non-negative. To show 
that this is the case, we note that the matrix 

N 

M[0} = l + l[B n (S13) 

n=l 

commutes with the following anti-unitary operator: 

U = is 2 a 3 K, (S14) 

where s are Pauli matrices which act on the spin index, a are Pauli matrices which act on the 
orbital (x,y) index, and K is the complex conjugation operator. Note that U 2 = —1. Using 
this, one can prove (26,27) (in a similar way to the proof of Kramers' theorem) that if \ a is an 
eigenvalue of M, A* is an eigenvalue also, and that if A a is real then it is doubly degenerate. 
The determinant can be written as det [M] = Y\ a l^af > 0. The integrand in the partition 
function (Eq. IS 121) is therefore non-negative, and can be simulated using Monte Carlo without 
a sign problem. Note that there are no particular restrictions on tij (e.g. it does not have to be 
bipartite) or fx. So particle-hole symmetry or specific densities are not required. 
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